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1 Introduction 

The Tokai to Kamioka (T2K) experiment uses a beam of muon neutrinos, produced at the 
J-PARC facility on the east coast of Japan, to study neutrino oscillations driven by the 
Ara^m mass splitting. A suite of near detectors located 280 nr from the secondary beam 
source samples the unoscillated beam, and the Super-Kamiokande water Cherenkov detector 
samples the beam at a baseline of 295 km, and at a point 2.5° off the beam axis, giving a 
narrow-band beam centred around 600 MeV. Analyses of the oscillation channels u^ —> u e 
and Vp —> u lt allow measurements to be made of 6 13 , $23 and Am^ tm , and, ultimately, for 
weak constraints to be placed on the CP-violating phase Scp ■ 

In addition to these analyses, T2K has made world-leading neutrino cross-section mea¬ 
surements in the sub-GeV energy range, utilising both the near and far detectors. The 
present work will discuss both the most recent measurements of the oscillation parameters, 
and these cross section analyses. 

2 Neutrino oscillations 

Observations of neutrino flavour change conclusively favour a model whereby the neutrinos 
have non-zero mass, and the mass eigenstates of the neutrino are not identical with the 
weak (flavour) eigenstates. Instead they are related by a unitary transformation, known as 
the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix C/pmns : 

M 

= GpMNS I u 2 I • (1) 

w 

Since a neutrino produced in a specific weak eigenstate u a is in a superposition of mass, 
and therefore energy, eigenstates, then after propagating some distance the different phases 
between the mass states will produce a superposition of flavour states. The neutrino may 
therefore be observed in a different flavour state up. The probability for a neutrino produced 
in a flavour eigenstate u a to be observed later in a flavour eigenstate up, depends on the 
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neutrino energy E. the propagation distance L, the mass splittings between the eigenstates, 
A rriy = m 2 — m 2 , as well as on the elements U a i of the mixing matrix, as 

P(y a ->• v p ) =5 a p 



This unitary transformation may be written as the product of three two-dimensional 
rotations with angles (# 12 , ^ 13 , # 23 ); with the addition of a complex phase 5cp , and CP- 
violation will be present unless 6cp is equal to 0 or 7 r. The additional degrees of freedom in 
the unitary transformation correspond to phases which are of no consequence for oscillation 
physics. Specifically, the matrix may be factorised thus: 

/10 0 \ / ci 3 0 si 3 e~ lScp \ / Ci 2 S 12 0 \ 

DpMNS = { 0 C 23 S 23 I I 0 1 0 j| — S 12 C 12 0 I , (3) 

\0 -S23 C 23 / \-Sl 3 e +i,5c,p 0 C 13 / \ 0 0 1/ 

where (s t j, c tJ ) correspond respectively to the sin and cosine of the angle dij. For the case of 
T2K, it is informative to write the oscillation probabilities in terms of these mixing angles, 
neglecting terms containing Amf 2 , since at T2K energies and baselines, Am^j -C 1. This 
leads to the simplified expressions 

p( En -> - 1 - sin 2 (2023) sin 2 (4) 

-t Ve) ~ sin 2 (^ 23 ) sin 2 (20i 3 ) sin 2 ■ ( 5 ) 

It can be seen that when considering only terms in Am^, the zy /t survival probability depends 
only on the angle 023, and the u e appearance probability on 0i3 and 622 ,. Dependencies on 
012 , 5cp, and extra terms due to the presence of matter, appear only at higher orders. 

3 The T2K Experiment 

T2K is a long-baseline neutrino oscillation experiment based in Japan pQ. A beam of muon 
neutrinos is produced at the J-PARC facility using 30 GeV protons from the Main Ring 
(MR) accelerator, in a conventional fashion. The protons are extracted and impacted onto 
a target to produce a secondary beam containing mostly 7r and K, which are then focused 
using a series of three magnetic horns, whose polarity can be changed to focus positive 
or negative hadrons, producing a dominantly or beam respectively. Neutrinos are 
produced in the decays of hadrons; the dominant process for a beam is (7r + —> p, + + 

Due to the presence of other decays and wrong-sign hadrons in the beam, an admixture of 
electron neutrinos, and wrong-sign muon neutrinos, is also present. 
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The T2K experiment is novel in that the main detectors are placed approximately 2.5° 
from the axis of the neutrino beam, giving a narrow-band beam with neutrino energies 
tightly focused around 600 MeV, the energy of the first oscillation maximum at the far 
detector. It also removes most of the high-energy tail of the beam, reducing the back¬ 
ground from high-energy deep inelastic and resonant n production events, which can be 
mis-reconstructed as quasi-elastic events at a lower energy. 

The far detector for T2K is the Super-Kamiokande water Cherenkov detector, at Mozumi 
Mine near Kamioka, on the west coast of Japan. Super-K is 295 km from J-PARC, and 
has a rock overburden of 1 km (2.7 km water equivalent), reducing non-beam backgrounds 
to a negligible level. It has a 50 kt total water mass, and 22.5 kt fiducial mass. The 
detector is cylindrical and split into an inner detector (ID) and outer veto region (OD), both 
instrumented with PMTs. Neutrino interactions are identified by observing the Cherenkov 
radiation from the resultant charged lepton; muons and electrons are distinguished by the 
different ring shapes produced, with muons producing a clean ring, and electrons a fuzzy 
ring due to showering. 

A near detector complex at the J-PARC site, 280 m from the target, is used to constrain 
the properties of the unoscillated neutrino beam. It consists of an on-axis component 
(INGRID) and a detector at the same off-axis angle as Super-K (ND280). The INGRID 
consists of thirteen lxlxl m 3 iron-plastic scintillator sampling calorimeter modules 
arranged in a plus configuration centred on the beam axis, with two additional off-diagonal 
modules. Its purpose is to measure the overall intensity and direction of the beam. 

The ND280 detector is a general-purpose detector comprising multiple components. At 
the upstream end there is a module (the POD) optimised for reconstruction of 7r° events. 
Downstream of this is a tracker region, consisting of two fine-grained scintillator detectors 
(FGDs) which make up the target mass, sandwiched between three gas TPCs. Surrounding 
the inner detectors are electromagnetic sampling calorimeters (ECALs), and a 0.2 T magnet 
to allow reconstruction of the momentum and sign of charged particles in the TPCs. The 
yoke of the magnet is instrumented with side muon range detectors (SMRDs), which also 
serve as a cosmic trigger. 

The results presented in this paper are based on data taken in the period 2010-2013 
(Runs 1-4). This period corresponds to a total of 6.6 x 10 20 30 GeV protons on target. 
T2K is presently gathering data for analysis. 

4 Flux constraint from beam simulation and ND280 

The nominal beam flux, including admixtures, is modelled using a Monte Carlo, which uses 
GEANT3 to track hadrons through the target and horns, and internally implements the 
relatively simple propagation and decay of particles in the beam pipe. Results of the NA61 
hadron production experiment at CERN are used to constrain hadron production in the 
target. 

Results from the ND280 are used to tune the nominal Monte Carlo. The analysis used 
for this tuning is based on identifying and classifying charged-current (CC) interactions 
of muon neutrinos in the most upstream FGD of the detector. z/^-CC events are tagged 
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by looking for a vertex in FGD1, where the highest-energy negative track in the TPC 
immediately downstream has energy deposition consistent with a muon. These events are 
then classified as CC-07T, CC-l7r + or CC-other, based on the pions identified in the final 
state. Charged pions are tagged by pion-like TPC tracks, or, for positive pions, Michel 
electrons in FGD1. Neutral pions are tagged using electron-like tracks in the TPC. Any 
event containing multiple pions, or a 7 r° or 7r _ , falls into the CC-other category. Figure [l] 
shows the momentum spectra for all events in these categories, both in data and nominal 
Monte Carlo. 



(a) v^CC-On (b) v^CC-ln (c) ^CC-other 

Figure 1: i/„-CC event spectra in the ND280, for data and nominal Monte Carlo 


The beam Monte Carlo parameters and cross-section models, with prior external con¬ 
straints, are fitted to the three ND280 samples, to calculate the expected event rate at 
Super-K. The resulting improvement in event rate predictions can be seen in Figure [2j 



Reconstructed v Energy (GeV) 

(a) 



Reconstructed v Energy (GeV) 


(b) 


Figure 2: Predicted event spectra at Super-K, before and after constraints from fitting to 
ND280 data. Shown for v e (a) and (b). 


5 Muon neutrino disappearance analysis 

Muon neutrino disappearance at T2K is observed via a deficit in the number of events 
seen at Super-K relative to that predicted in the absence of oscillations [2|. This deficit is 
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expected to be energy-dependent, and enables measurements of 023 and Am| 3 (Am| 1 ) for 
normal(inverted) hierarchy (i.e. the sign of the mass splitting Am^), hereafter referred to 
as NH or IH. 

A cuts-based analysis is used to obtain an event sample enriched with CCQE events 
{yn + n —> + p). The proton is usually below the Cherenkov threshold and therefore 
invisible. The analysis demands a single ring, with muon-like PID, vertex and whole event 
in the fiducial volume, low activity in the OD, visible energy > 30 MeV, reconstructed 

> 200 MeV, and < 1 decay electron. The event must be in time with the beam. For 
selected events, the neutrino energy is reconstructed assuming that the underlying neutrino 
interaction is quasi-elastic. 

Oscillation parameters are estimated from this selection using an unbinned maximum- 
likelihood fit to the events. A full three-flavour fit, including matter effects, is performed. 
The oscillation parameters Am^Am^) for NH(IH), Scp and sin 2 623 are allowed to float 
freely in the fit. sin 2 $ 13 , sin 2 $12 and Am^ are fitted with external constraints from reactor 
and solar+KamLAND data. 45 systematic parameters are included in the fit, from cross 
section uncertainties common with ND280 and specific to Super-K, Super-K detector ef¬ 
fects, final state (FSI) and secondary (SI) interactions, and the aforementioned constrained 
oscillation parameters. The systematic errors from these sources are broken down in Table 

m 


Source of uncertainty (number of parameters) 

c exp / exp 
5 n SK/ n SK 

ND280-independent cross section (11) 

4.9% 

Flux & ND280-common cross section (23) 

2.7% 

SK detector & FSI+SI systematics (7) 

5.6% 

sin 2 0 i 3 , sin 2 0 i 2 , Am^, 8 C p (4) 

0 . 2 % 

Total (45) 

8 . 1 % 


Table 1: Effect of la systematic parameter variation on the number of 1-ring /r-like events, 
computed for oscillations with sin 2 023 = 0.500 and |Am| 2 | = 2.40 x 10 ” 3 eV 2 /c 4 . 

Figure [ 3 ] shows the central values and 2D confidence regions for 023 and Am| 3 (Am| 1 ) for 
NH(IH), along with 90% C.L. limits for MINOS and Super-K for comparison. Confidence 
regions are calculated using the Feldnran-Cousins technique. The best fit values and ID 
68% C.L. intervals are sin 2 023 = 0.5 14^QQ 3 g(0.511 ± 0.055) and Am 2 3 = 2.51 ± 0.10 x 
10” 3 eV 2 /c 4 (Am 31 = 2.48T0.10 x 10” 3 eV 2 /c 4 ) for NH(IH). The error budget is dominated 
by statistical errors. 

6 Electron neutrino appearance analysis 

To measure electron neutrino appearance, the number of electron neutrino events, and their 
spectral shape, are compared to those expected in the absence of oscillations [3]. Background 
events come from the intrinsic v e component of the beam, and also from mis-reconstructed 
vp events, in particular NC- 7 T 0 + p —> + p + 7 r°), where one 7 from the pion decay is 
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Figure 3: Oscillation best fit point and confidence regions for the T2K anal¬ 
ysis, with MINOS and Super-K results for comparison. External constraints on 
other oscillation parameters are: sin 2 0 i 3 = 0.0251 ± 0.0035, sin 2 0 i 2 = 0.312 ± 0.016, 
Am^i = (7.50 ± 0.20) x 10 ~ 5 eV 2 /c 4 . 5cp is unconstrained and a uniform Earth matter 
density of 2.6 g/cm 3 is used. 


missed. The analysis allows measurement of sin 2 2013. 

The u e selection at Super-K is made using similar preselection cuts as for the v ]t analysis, 
on event timing, vertex location and event containment, visible energy and OD activity. As 
with the v,, analysis, a CCQE-enriched sample is obtained by selecting events containing 
a single electron-like ring, with reconstructed momentum greater than 100 MeV/c. Events 
with a reconstructed neutrino energy (assuming CCQE kinematics) greater than 1250 MeV 
are also rejected, since these are dominantly due to the intrinsic u e beam component. To 
remove 7r° background, an algorithm measuring the maximum likelihood of the observed 
PMT hit pattern, assuming zvCCQE ca NC- 7 T 0 , is used; events are accepted or rejected 
based on the reconstructed 7r° mass and the likelihood ratio between the 7r° and electron 
hypothesis best fits. These cuts lead to an expected 21.6 selected events for sin 2 2013 = 0.1 
and 5cp = 0; the expected number for $13 = 0 is 4.92. In actuality, 28 candidate events 
were observed. The total systematic error on the number of observed events is 8 . 8 % for 
sin 2 2013 = 0 . 1 ;this is dominated by neutrino cross-section uncertainties. 

A binned extended maximum-likelihood fit is used to extract oscillation parameters 
from the data. The likelihood function contains terms for the signal shape, signal nor¬ 
malisation, T2K systematics and external oscillation constraints. Events are parameterised 
using the reconstructed electron momentum and angle; however, parameterising events 
using reconstructed neutrino energy gives compatible best-fit points and near-identical 
confidence regions. Other oscillation parameter inputs are as follows: sin 2 dio = 0.306, 
Am. 2 i = 7.6 x 10 ~ 5 eV 2 , sin 2 023 = 0.5, |Am 2 2 | = 2.4 x 10 ~ 3 eV 2 , 5cp = 0. For the NH(IH) 
the best-fit values, and 68 % confidence level bounds, are sin 2 2013 = 9.140^q ^(O.^O^QQg 3 ). 
The significance of non-zero 0 i 3 is found to be 7.3er; the same result is obtained with ei¬ 
ther a delta log-likelihod method or toy Monte Carlo. This significance was obtained using 
the values for sin 2 023 and 5cp as above; however the significance remains above 7.0cr for 
any values of these parameters consistent with their uncertainties. The sin 2 2013 confidence 
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limits are shown as a function of 5cp in Figure |4j 
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Figure 4: Oscillation best fit point and confidence regions for the T2K u e analysis. 


7 Cross section measurements 

The T2K experiment is making world leading neutrino cross section measurements in the 
~ 1 GeV range. Only results and outlines of the analyses will be presented here; the reader 
is directed to the references given for detailed descriptions. 

7.1 u e CC-inclusive cross section on Carbon 

The CC-inclusive cross section of v e on Carbon has been measured in the ND280 detector [3] - 
The signal for this process is any event containing an electron in the final state, and the 
principal background is 7 s from 7T° decay. 

Events are selected which have a negative electron-like track (by dE/dx) in TPC2, with 
a vertex in FGD1. Events are subject to a veto on activity upstream of FGD1. The 7r° 
background is reduced by demanding that events containing a positive track do not have an 
invariant mass compatible with a 7 r° when the track is matched with the electron candidate. 

Figure [ 5 ] (a-c) show the differential cross sections (averaged over the T2K flux), as a 
function of electron momentum, interaction Q 2 , and the cosine of the electron polar angle, 
respectively. Good agreement is seen between T2K and the generators for the differential 
cross sections. Figure [5] (d) shows the measured cross section, integrated over the T2K 
beam flux, along with results from the Gargamelle bubble chamber experiment and the 
NEUT and GENIE Monte Carlo interaction generators; again, good agreement is seen with 
the generators, and also the Gargamelle data. 
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Figure 5: z/ e -CC cross sections on Carbon, integrated over the T2K beam flux, (a-c) show 
differential cross sections in electron momentum, interaction Q 2 , and the cosine of the 
electron polar angle, respectively, (d) shows the total cross section. 


7.2 CC-inclusive cross section on Carbon 

The ND280 detector has also measured the z/ /t CC-inclusive cross section on Carbon [5]. 
The selection for this measurement is essentially identical to that used for the analysis con¬ 
straining systematics for Super-K (see Section [4]). All z/^-CC candidate events are included 
in the selection. 

The events are binned in muon momentum p^ and polar angle cos 8 ^, and differential 
cross sections are computed in these values, integrated over the beam flux. These results 
are shown in Figure [6]— the results agree well with the NEUT and GENIE generators. 

7.3 NCQE cross section from de-excitation gammas 

T2K has made the first ever measurement of neutral current quasi-elastic (NCQE) neutrino 
scattering on Oxygen [bj. This process, v x + 16 0 —> u x + p + 15 N* (or u x +n + 15 0*), leaves 
the daughter nucleus in an excited state, so the events can be tagged using gamma rays 
from nuclear de-excitation. A beam with a low duty cycle, combined with precise timing 
information, allows to reduce the non-beam background to a negligible level, enabling a 
measurement which would be impossible with atmospheric neutrinos. 

After making standard Super-K cuts to remove non-beam background, and vetoing on 
high-energy preceding events, events are selected which have a visible energy between 4 and 
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Figure 6 : z/^-CC cross sections on Carbon, integrated over the T2K beam flux, (a-d) each 
show the differential cross section in muon momentum, for a different range of muon polar 
angle. 


20 MeV, and a large Cherenkov angle. This Cherenkov angle cut helps to remove low-energy 
muons from the event sample. After making these cuts, the dominant background is from 
NC non-QE events. 

Figure [7] shows the measured cross section, averaged over the T2K beam flux. This is 
seen to be compatible with the Ankowski NCQE model over the flux energy range. 

8 Summary 

The T2K experiment has published results showing exclusion of $13 = 0 of 7.3 a in the v e 
appearance channel, and giving the world’s most precise measurement of # 23 - It has also 
produced world-leading neutrino cross section results using both near and far detectors. 
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Figure 7: Measurement of the NCQE cross section for the T2K beam flux. 
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